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Abstract. We study the dynamics of billiard models with a modified collision rule: the 
outgoing angle from a collision is a uniform contraction, by a factor A, of the incident 
angle. These pinball billiards interpolate between a one-dimensional map when A = and 
the classical Hamiltonian case of elastic collisions when A = 1 . For all A < 1 , the dynamics is 
dissipative, and thus gives rise to attractors, which may be periodic or chaotic. Motivated 
by recent rigorous results of Markarian, Pujals and Sambarino |MPS08|, we numerically 
investigate and characterise the bifurcations of the resulting attractors as the contraction 
parameter is varied. Some billiards exhibit only periodic attractors, some only chaotic 
attractors, and others have coexistence of the two types. 



AMS classification scheme numbers: 37D45, 37D50, 37M25, 65P20 

1. Introduction 

Billiard models have played an important role in mathematical physics for forty years, since 
their introduction by Sinai in [SinTOJ, motivated by the Boltzmann ergodic hypothesis on 
fluids of hard spheres. They are a class of models which are accessible and of interest 
both to mathematicians |CM061 and to physicists IIGas98i In particular, they have attracted 
much attention as simple models of physical systems which exhibit strong chaotic properties 
(hyperbolicity, ergodicity and mixing) 

In classical billiards, particles travel freely between collisions with hard, fixed 
boundaries. They then collide elastically, i.e. without friction, so that the kinetic energy 
is conserved. Such billiards can thus be considered as a type of Hamiltonian system with 
discontinuities. 

Our objective is to investigate, mainly numerically, the different types of dynamics which 
occur in billiard tables with a modified collision rule: we consider pinball billiards, where the 
reflection rule at collisions is no longer elastic. Rather, the angle (with respect to the normal) 
of reflection is smaller than the angle of incidence, by a constant factor The dynamics is 
then no longer conservative; rather, there is contraction in phase space, which gives rise to 
attractors. Two types of attractors are found to occur: they are either (i) periodic points whose 
stability matrix has eigenvalues with modulus not greater than one, or (ii) non-trivial invariant 
sets that exhibit dominated splitting (defined below). Both types of attractors can coexist, and 
we study examples where all possible cases occur. 

Rather than attempt to model a particular physical scenario, we consider a model that 
interpolates between a one-dimensional dynamics, where reflection is always normal to the 
boundary, and a Hamiltonian one (elastic reflection). A related modified billiard collision 
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rule, which was motivated by physical considerations, was studied recently in IIADMH08I . 
Other Hamiltonian systems with dissipation have been studied from a similar point of view, 
for example the kicked rotor map IIFGHY96I IMG08I . 

We are interested in the dynamical properties of orbits in the forward limit set of the 
system, that is, the union of all OJ-limit sets of points in the space, which describes the 
dynamics of the system as the time goes to +00. Markarian, Pujals and Sambarino IIMPS08II 
proved that, under certain circumstances, case (ii), that is the ocurrence of non-trivial invariant 
sets with dominated splitting, appears in many common billiard tables. In our context, an 
invariant subset A of the limit set has dominated splitting if the tangent bundle continuously 
splits into two non-trivial invariant sub-bundles under the derivative, say T/^ — E(BF, and 
there are C > and v G (0, 1 ) such that 

ll^rifiM II • ||i).r"if (/'.(.)) IK cv" (1) 

for any x €E A and n ^ 0. In these circumstances one can expect fMPSOSl that the limit set can 
be decomposed into the union of three subsets: (1) a finite union of disjoint invariant curves 
formed by periodic points with finite period; (2) a finite union of simple closed curves where 
the dynamic is conjugate to an irrational rotation; and (3) a disjoint union of finitely many 
compact invariant hyperbolic sets. 

In this paper we restrict attention to tables with focusing boundaries. A focusing 
boundary (which has negative curvature under our convention) has the property that a beam 
of parallel trajectories which hits it focuses after the collision. Bunimovich showed I Bun79ll 
that if there is "enough space" in the billiard table, then after it has focused, the beam can then 
defocus, before it collides again. This defocusing mechanism can give rise to strong chaotic 
properties, in particular hyperbolicity |Woj86 IMar88i and to systems which are ergodic and 
mixing, see for instance IISza92l and | MarQSf! 

One of the motivations of this work is to study pinball billiard maps defined on tables 
which in some sense are close to the limit of validity of the hypotheses assumed in |,MPS08l 
for which rigorous results could be obtained. We consider cases in which either the attractor 
is not contained in a compact region of the phase space, or certain geometrical properties used 
in that reference do not hold. In this way, we shed light on the interesting types of behaviour 
which pinball billiards can exhibit. 

We begin by defining formally the dynamical system, and then proceed to study 
numerically the dynamics in several different types of billiard table with focusing (and in 
one case neutral) boundaries. The aim of the paper is to characterise the types of attractor that 
can be expected to occur in pinball billiards, and relate them to the geometric features of the 
billiard table. In particular, there is a competition between the defocusing mechanism and the 
angle contraction in the pinball dynamics. 

The tables we consider are the family of ellipses, the cardioid, a modified, convex version 
of the cardioid with a discontinuity in the curvature, and a smooth family of perturbations of 
the circle. In the ellipse, it is found that the limit set is always periodic and it can either have 
dominated splitting or not, depending on whether certain eigenvalues are real or complex. 
In the cardioid, the limit set always consists of a chaotic attractor which accumulates on the 
vertex. We then introduce a modified, convex version of the cardioid, in which the vertex 
is removed. There we find that the limit set consists of periodic or chaotic attractors, and 
for some values of the contraction parameter these two phenomena coexist. Here there are 
two points where the curvature is discontinuous, and the chaotic attractor contains trajectories 
which include these points. Finally, in a smooth family of perturbations of the circle we 
observe coexistence of periodic and chaotic attractors, which undergo bifurcations when the 
contraction parameter is varied. 
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2. Pinball billiards 

Let Q be an open, bounded, simply-connected subset of the plane, whose boundary Y consists 
of a finite number of compact curves F,. Billiard dynamics in Q is the dynamical system 
describing the free motion of a point mass particle inside Q with reflections at the boundary 
r :— U,r/. Each F, is called a (smooth) component of F. Let n{q) be the unit normal vector 
at the point q E F, which points towards the interior of Q. The configuration space of this 
dynamical system is Q, and its phase space M is given by 

M:={{q,v):qE F, |v| = 1, (v,n(<?)) ^ 0}, (2) 

where (•) denotes the Euclidean inner product. 

We introduce coordinates {s, (p) on the phase space M, where s is the arc length parameter 
along F and (p is the angle between the reflected vector v and the normal n{q), satisfying 
(v, n{q)) = cos (p. In classical billiards the angle of incidence is equal to the angle of reflection 
(elastic reflection) |,CM06]; in this case, the billiard map is area-preserving in terms of the 
coordinates (i,sin((j!>)), known as Birkhoff coordinates. 

Recently, Markarian, Pujals and Sambarino IIMPS08II introduced pinball billiards, which 
are billiard maps with a modified reflection rule. Let rj E [— f , f ] be the angle from the 
incidence vector to the outward normal —n{q) at the point q where the ball hits the boundary. 
In pinball billiards we no longer restrict to the elastic case; rather, we allow the outgoing angle 
(p to depend on the incidence angle 7] . In principle, a dependence could be allowed on the 
position of the conision, as in BADMH08I . but we suppose that the exit angle (p depends only 
on the incidence angle T]. Then (jd — —f]+f{f]), where / : [—n/2,n/2] ^ M is a function. 

In [MPS08J, section 3.2 of which gives more details on these topics, two perturbations 
on the angle of reflection were studied. In this paper, we restrict attention to one of them: 
we assume that f{r\) — (1 — A))], with constant A < 1; thus <p — — At], that is, the outgoing 
angle is given by a uniform contraction of the angle of incidence. For simplicity we rescale 
the speed of the particle after each collision so that it is always I. Physically this corresponds 
to a sticky boundary which gives exactly the right amount of "kick" to the particle; this can 
be thought of as a caricature of what happens in a real pinball table, whence the term "pinball 
billiard". 

The pinball billiard map Tx with a given contraction parameter A E [0, 1] is then defined 
as follows: we leave the point G T with direction vq until we intersect the boundary again 
at the point qi E F. If the unit normal n{qi) at the new position is well defined, then the new 
velocity vector vi is given by the above collision rule. The pinbilliard map is then defined 
by Tx{qo, Vq) :— {qi,vi). lfn{qi) is not uniquely defined, for example if the trajectory lands at 
a sharp comer, then the map is not defined at (^o, t^'o)- ?! is also undefined at points {qo,vo) 
where the trajectory hits tangentially another component of the boundary F. Removing such 
points from the domain of definition of the map forces us to work in a non-compact domain. 
This is an obstacle for the direct application of results from BMPSOBI . 

Let f be the Euclidean distance between qo and qi , and Ki be the curvature of F at qi taken 
with the following convention: Ki is negative if the component of the boundary is focusing, 
zero if it is flat, and positive if it is dispersing. Zero curvature is not allowed on focusing 
and dispersing components of the boundary. In this paper we will only study billiards with 
focusing and flat components, and only for A ^ 1. Although it is possible to consider A > 1, 
in this case the trajectories tend to converge to "whispering gallery" modes tangent to the 
billiard boundary. 
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2.1. Derivative matrix 

The total derivative of the pinball billiard map Tx at a point x — (5, (po) G M is, as shown in 
BMPSOSI subsection 3.3], given by: 

D^Us^cp,) = - ^ ^K,l+K,) X {K^B+l) ) ' 

where 

tKo + cos (Dn t 
COSTJi COST]! 

Again, subscript refers to quantities evaluated at the starting collision point, and subscript 1 
to quantities at the next collision. Thus <Pq is the exit angle of the trajectory at the boundary 
point from which the trajectory leaves, and 771 is the incidence angle at the next collision, 
calculated by the rule described above. The derivative thus involves information from two 
consecutive collisions. 

The stability matrix of an orbit is the product of the stability matrices of each trajectory 
segment. In the case of a periodic orbit, it is sufficient to take the product over a single period. 

2.2. Slap billiard map: A = 

The pinball billiard map depends on the contraction parameter A, and one of our goals is to 
study the global modifications (bifurcations) of the dynamics of Tx as A varies in the interval 
[0, 1]. When A = 1 the outgoing angle is equal to the incidence angle, so that T\ is the classical 
elastic billiard map, and the measure djU = cos ^d.sd^ is preserved. 

When A = 0, trajectories leave exactly along the normal vector at each collision, so that 
Tq reduces to a one-dimensional map, which was named the slap billiard map in IIMPS08I : 
the derivative of 7b at a point q along the boundary is given by T^^q) = — (fA'o + l)/cosT7i. 
Note that the slap billiard map is in fact a map of the circle, since the arc length is a periodic 
variable. We emphasise that in the cases considered in this paper the curvatures Ki are always 
non-positive. 

We remark that the slap billiard map is related to the geometrical concept of evolute of 
a curve, which is the envelope of the set of all normals to the curve. Since the trajectories in 
the slap billiard always follow the normal vectors, the evolute is generated in configuration 
space. Nonetheless, this seems not to be particularly relevant for the dynamical properties of 
the map. 

3. Circular and elliptical tables 

The simplest smooth, convex tables with focusing boundaries are elliptical tables, which are 
classical examples in the case of elastic billiards IICM06I Chapter 1]. In this section we study 
pinball billiards in elliptical tables, beginning with the special case of circular tables. The 
tables have elliptical boundaries, parametrised as Ta{t) — (acosf ,sinf), where t G [0,2;r) and 
the shape parameter a satisfies a ^ 1 without loss of generality. 

3.1. Circular table 

The simplest case of pinball billiard dynamics occurs in a circular table, for a = 1 . When 
A = 1, that is, for elastic collisions, the incident and outgoing angles are constant along 
each trajectory, due to conservation of angular momentum, and the phase space is foliated 
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Figure 1. Single trajectories in configuration space for pinball billiard dynamics with 
contraction parameter A = 0.99 in (a) circular table, and (b) elliptical table with a = 1.5. 
Colours indicate the number of bounces, with lighter colours corresponding to later times, 
exhibiting the asymptotic convergence to period-2 orbits. The initial condition in (a) is a 
random one; that in (b) was taken close to the unstable period-2 orbit along the major axis, 
from which it rapidly diverges. Remnants of the structure of caustics which is found for the 
integrable case (X = 1) are visible in (b). 



by parallel lines corresponding to constant angles. This system is thus completely regular 
(integrable). 

For A < 1, however, at each collision the exit angle (p' = X(p decreases, so that the 
outgoing trajectory is closer to the diameter which corresponds to the normal direction at 
the point of collision. This continues at each collision; since there is no counter-effect, the 
trajectory from any initial condition converges to some period-2 diametrical orbit, as shown 
in figure [TJ a). Any diameter is a limit set, and the exact limiting behaviour depends on the 
initial conditions. 

In IIMPS08I it is proved that for < A < 1 the limit set of the pinball billiard in a circular 
table has dominated splitting. This limit set consists of the line [(p = 0], the set of all period-2 
orbits. 

The slap billiard map of a circular table corresponds to a translation by n on the circle: 
To{s) = s + 7Z (mod In), that is, each point s is mapped to the diametrically opposed point 
s + n and then back, so that each point is periodic with period 2. 

3.2. Elliptical table 

New dynamical features appear when we break the circular symmetry and consider an ellipse 
with non-zero eccentricity. An ellipse has exactly two periodic orbits of period 2, lying along 
its major and minor axes. However, these two periodic orbits are found to have different 
stability properties under the pinbilliard dynamics: the orbit lying along the minor axis is 
stable, whereas the orbit lying along the major axis is unstable: an initial condition pointing 
along the major axis stays there forever, whereas an initial condition starting arbitrarily close 
to it leaves the neighbourhood of the major axis and eventually converges to the stable orbit, 
as shown in figure[TJb). 

The complete description of the dynamics is that there is one stable and one unstable 
period-2 orbit and they form the limit set of the pinball billiard transformation, for all 
A e (0, 1). With the exception of the unstable orbit, all other trajectories are attracted by 
the stable one. 
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(6) 



3.3. Stability of period-2 orbits 

To confirm the results on the stability of the period-2 orbits, we calculate explicitly their 
stabihty matrices. The curvature of the ellipse boundary at the point r„(r) is given by 

K{t) - (5) 

[a2sin2(f)+cos2(f)] ' 

Thus the curvature at the points of the ellipse lying on the major axis is K'(O) = —a, and on 
the minor axis is K{n/2) = —l/a^. 

The stability matrix of the periodic orbit of period 2 along the minor axis, for contraction 
parameter A, is the product of the stability matrices for the two flights of the orbit, and is 
found to be 

/ [4(l+A)(l-fl2)^_^4j 2(l+A)(fl2_2)/fl2 \ 

V -2A(l+A)(fl2-i)(fl2_2)/fl6 A[4(l-fl2)(i+A)+Afl4]/fl4 j' 

The stability of the orbit is determined by the eigenvalues of (j6|l, which are given by 

Pl{a,X)±^/PA^) 
M±(a,A) = , (7) 

where 

Pi = - 2)2A2 - 8(fl2 - 1)A + (fl2 _ 2)2; 

P2 = [(1 + A ) (a^ _ 2)]2 [(fl2 _ 2)2a2 -2{a* + Aa^ - 4)A + {a^ - 2)2] . (8) 

The eigenvalues are complex conjugates satisfying (A ) | < 1 for A close to 1 and any a^ I, 
so that the minor axis is then a stable focus. Furthermore, |/z±(A)| 1 when A ^ 1, since 
the orbit is elliptic in the elastic case A = 1 . 

The eigenvalues change type from complex to real when P2{a,X) crosses 0. This occurs 
when a = \/2, for any A e [0, 1]. For a ^ \fl, we have also the following solutions: 



a 



' 4fl2-4±4fl2Vfl2-l 



A±(«) = . (9) 

We find that A+ is always outside the region of A of interest, but A_(fl) £ [0, 1] for all a ^ 1, 
with A„ (a) when a \/2 and A_ (a) 1 when a — * oo. For pairs (a, A ) below the graph 
of A_(a), both /z±(fl,A) are real and the dynamics is then that of a stable node; above this 
graph, both eigenvalues have non-trivial imaginary part. These eigenvalues are shown for the 



case fl = 5 in figure 2(a) Note that the absolute value of the eigenvalues must reach 1 when 
A = 1, since the orbit is elliptic in the elastic case. 

The stability of the other period-2 orbit, along the major axis, can be found by inverting 
the geometry. A computation shows that both eigenvalues jj.± are real, with < jU_ < 1 < 
Hence the major axis is always an unstable period-2 orbit for A < 1. 

In our two-dimensional setting, the first requirement to have a set with dominated 
splitting is the existence of two different non-trivial invariant directions for Z)7\. Thus, if 
the limit set contains a periodic point with non-real eigenvalues, then this set cannot have 
dominated splitting. On the other hand, if the limit set consists of only a finite number 
of hyperbolic periodic orbits with real eigenvalues, then the limit set does have dominated 
splitting. Hence, for any a ^ 1 there is a number Xa such that the limit set of has dominated 
splitting if A e (0, Aa] and does not have dominated splitting for A G (A„, 1). 
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3.4. Slap billiard map of an elliptical table 



The slap billiard map To for an elliptical table with a = 1.2 is shown in figure 2(b) together 
with its second iterate. We use the polar angle Q of the collision point instead of the arc length 
along the ellipse, since it is simpler to calculate. In the configuration space (not shown), 
all trajectories of the slap billiard map are tangent to the evolute of the ellipse, which is a 
stretched astroid (the Lame curve). 

On the graph of Tq can be seen the two unstable period-2 points lying on the major axis 
(0 = and Q — n), and the two stable period-2 orbits on the minor axis. The slap billiard 
map of any elliptic table is not an expanding map, which agrees with the absence of chaotic 
behaviour in the dynamics. 



4. Cardioid 

The results of IIMPS08I suggest that pinball billiards may have much more complicated limit 
sets than just attracting periodic orbits, in particular non-trivial invariant sets with hyperbolic 
dynamics. From here on we search for more complicated dynamics. 

In this section we study pinball billiards in the cardioid, parametrised by the equation 
p{0) = 1 +cos(0) in polar coordinates (p , 0), with —n^9<n. The standard, elastic billiard 
dynamics in a cardioid has been studied extensively |Woj86 |Mar93l ISza921 'BD97 DBOll 



and in particular has been proved to be hyperbolic, ergodic and mixing. Robnik |Robn83l 
studied the existence, stability and bifurcations of periodic orbits in a class of billiards which 
includes the cardioid, with the focus being the transition from regular dynamics in the circle 
to stochastic dynamics in the cardioid. Here we rather study the dynamics in this single 
table when we vary the contraction parameter A between the one-dimensional map and area- 
preserving cases. 

For the cardioid, the arc length s{6) up to an angle 9, measured from the positive 
jc-axis, can be evaluated exactly as s{9) = 4sin(0/2) for ^ 9 ^ n. Thus i = is the 
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rightmost (flattest) point of the cardioid, positive values of s correspond to the top half of 
the cardioid, and the cusp is at s — ±4. This s is used as the abscissa in the phase space plots 
below. Some details of the numerical algorithm used to follow the dynamics are given in the 
appendix. The curvature of the boundary may be evaluated using the standard expression in 
polar coordinates: 

yM±im^fW^. (10) 
\p{<Sf + P'{»)-\' 

The negative sign here corresponds to our convention that focusing boundaries have fc < 0. 
The curvature of the cardioid is then K{d) ~ — 3/-\/8(l +cos 9). 

4.1. Chaotic attractor 

Numerical evidence shows the existence of a non-trivial chaotic attractor for any X E [0, 1), 
as shown in figure [3] The attractor grows in phase space as X increases, until it fills the whole 
of phase space uniformly when A = 1 . Here and in the rest of the paper, the angular coordinate 
in phase space is shown as sm{(p). 

We remark that points which are arbitrarily close to the singularity = ±7r at the cusp 
of the cardioid (corresponding to s ~ ±4) seem to belong to the attractor However, since 
the cusp does not belong to the domain of the pinball billiard map, since the normal vector is 
not uniquely defined there, the attractor does not lie in a compact region of phase space. 
Thus Theorem 5 of |MPS08| cannot be automatically applied. However, the dominated 
splitting behaviour foreseen by Theorem 4 of that reference is suggested by the following 
computations. 

The trajectory emanating from (almost) any initial condition is found numerically to have 
two Lyapunov exponents v±, with v+ > > V_, i.e. they are uniformly bounded away from 
zero for any A € (0, 1]. The numerical calculation of these Lyapunov exponents is summarised 



in the appendix, and the results are shown in figure 4(a) Due to the dissipative nature of 
the dynamics, there is contraction in phase space, which implies that v+ + V_ < 0. Each 
point on the graph, for a fixed value of A, is an average over 100 initial conditions. The 
Lyapunov exponents evaluated along the different trajectories are found to vary little - error 
bars are of the size of the symbols. This, together with visualisations of the attractor for 
different initial conditions, suggests that there is transitive behaviour on the attractor, i.e. that 
any initial condition on the attractor visits the whole attractor. These numerical results lead 
us to conjecture that the system is in fact non- uniformly hyperbolic in the sense of Pesin, 
i.e. almost every point has non-zero Lyapunov exponents IIRobi99l Section 9.2]. However, 
rigorous proofs of such conjectures are notoriously difficult IIRobi99L 

4.2. Slap billiard map of the cardioid 

In configuration space, all trajectories of 7b are tangent to the evolute of the cardioid, which 
is a mirror image of the cardioid scaled by a factor of 1/3 (not shown). The slap billiard 
map of the cardioid is interesting in that it has regions, close to the point i = of minimal 
curvature, where the derivative Tq is less than one, and other parts where it is greater than 



one, as shown in figure 4(b) Nonetheless, as the figure shows, the second iterate Tq of 
this map is piecewise expanding, i.e. its derivative is uniformly bounded away from 1, with 
(Tf^yis) > 1 everywhere. This confirms the strong chaotic properties of the dynamics, and 
suggests a situation analogous to that in Theorem 12 of 0MPSO8I with focusing walls, which 
states that under these conditions there should be a finite number of expansive attractors. In 
fact, numerically we find that the dynamics of the map seems to be ergodic. 
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Figure 3. The chaotic attractor in phase space for the cardioid, for (a) A = 0.3 and (b) A = 0.8. 
The inset of (a) shows the attractor in configuration space. The coordinates are arc length s 
and sin((p), where <p is the exit angle at each collision. The cusp of the cardioid is at s = ±4. 




Figure 4. (a) The two Lyapunov exponents v± for the cardioid, as a function of X . Their sum, 
which is the logarithm of the area contraction in phase space under the map, is also shown. It 
converges to — oo when A — > 0, since the map reduces to a one-dimensional map in this case, 
and to when A ^ 1, the elastic case, (b) The slap billiard map Tq for the cardioid and its 
second iterate. To itself has regions where the derivative Tq is less than 1 in modulus, near the 
cusp at s = 0. Nonetheless, its second iterate is expanding, with |rQ'(i)| > 1 everywhere, and 
thus the dynamics of the map are strongly chaotic. The two unstable fixed points of the second 
iterate correspond to an unstable vertical period-2 orbit joining the top and bottom points of 
the cardioid at 6 = ±;r/3, i.e. at i = ±2. 
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5. Cuspless cardioid 

In the previous two sections, we have studied tables which show completely regular and 
completely chaotic behaviour, both for elastic and pinball billiard dynamics. In this section 
and the following one, we consider two convex tables which exhibit a complicated mixture 
of phenomena. Here we begin with a convexified version of the cardioid, obtained by 
flattening the cusp of the cardioid. To do so, we remove the part of the cardioid with 
d e [—n, — ^] U and replace it by a vertical Une at x = — ^ between y — ±^\/3. 

This procedure gives a convex table which is smooth at all points except for these two joins, 
where the curve is only and has discontinuous curvature. 

Again the arc length can be evaluated exactly. We now take = as corresponding to 
the centre of the vertical line at the left of the table, and the rightmost point being at 5 = ±C, 
where C := 4sin(;r/3) + V3/4 = 9^3/4. The discontinuities are then at s = ±%/3/4. 

By unfolding this "cuspless cardioid" by reflection in the vertical line, we obtain an 
equivalent table which satisfies Wojtkowski's defocusing condition |Woj86 IMar88l . and 



is thus hyperbolic and ergodic in the classical case of elastic collisions iSza92l IMar93l . 
However, this unfolding process is no longer valid in the case of pinball billiards, since the 
outgoing angle is modified at each collision with the vertical line. 

From the point of view of pinball dynamics, the dynamics always has dominated splitting 
IIMPS08[ Theorem 4]. This is verified later, by computing Lyapunov exponents. However, 
once again the phase space is not compact. Thus the results of LPS 09 J can be immediately 
applied only if the attractor is contained in a compact subset of the phase space. This is the 
case for X close to zero. However, a non-trivial attractor is observed only for larger values of 
X. In fact, it is related to the existence of points of the boundary which are not C^: generic 
orbits in the chaotic attractor seem to accumulate on these non-smooth points. 

Numerical experiments show that for A close to 0, the period-2 orbit lying in y = 
attracts almost all initial conditions, except for the unstable vertical period-2 orbit which is the 
same as in the unmodified cardioid. For X close to 1, however, we observe a globally attracting 
chaotic attractor, as in the unmodified cardioid, although rather different in structure from that 
case. Careful simulations show that in fact there is a region where the period-2 and chaotic 
attractors coexist: some initial conditions are attracted to the periodic orbit, and others to the 
chaotic attractor. Figures |5] and |6] show the evolution of the numerically-determined attractors 
as X is varied, in configuration space and in phase space, respectively. In the following we 
shall explain these phenomena. 

5.7. Period-2 orbit 



We start by studying the horizontal period-2 orbit which runs along the symmetry axis y — 0, 
shown in figure |5j a). This orbit exists for any X £ [0,1], but its stability depends on the 
eigenvalues of the stability matrix, i.e. the derivative of at the point x = (0,0), which in 
turn depend on X. The stability matrix is obtained by composing the stability matrices at the 
points X and Ti{x) — (C,0), where C is the arc length to the point opposite 0. We find that 

„ 7^2_1/^ -4(27A + 11) 144(A + 1) \ 

'^■"'^-My 33X{x + i) -4X{nx+2i) )■ 

The eigenvalues ji of this matrix are 

H±{X) = ^ (^-11 -54A - 11A2±^(1 + A)v/11 + 86A + 11A2^ .(12) 

Both eigenvalues are real and satisfy jU±(A) ^ 0, for any A G [0, 1]. Moreover, ji±{X) > —I for 

'si 



any A e (0,Ac), where Ac := 5(27-8^11) ^0.0934 is the value such that that jU-(Ae) = -1, 
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(a) A = 0.02 



(b) A = 0.072 



(c) A = 0.2 



(d) A = 0.5 



1.0 



0.5 



0.0 



-0.5 



-1.0 



Figure 5. Evolution of numerically-observed attractors in configuration space with increasing 
A. ForA < A, ~ 0.0712, there is just aperiod-2 attractor. This periodic attractor coexists with 
a chaotic attractor for A e [A,,Ac], where A^. ~ 0.093. The period-2 attractor then becomes 
unstable, leaving just the chaotic attractor, which expands for increasing A . 
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(a) A e [0.1,0.3] 



(b) A e [0.4,0.6] 



(c) A e [0.7,0.9] 



Figure 6. Evolution of the chaotic attractor of the cuspless cardioid in phase space, as a 
function of A. The attractor is depicted for A between 0.1 and 0.9 in intervals of 0.1, with 
three values of A shown in different colours in each figure (red, blue and cyan in increasing 
order of A). Vertical lines mark the centre of the vertical section of the boundary at i = and 
the two curvature discontinuities at i = ±1/3/4. The attractors are also seen to be non-smooth 
here; they expand in phase space as A increases. 



i.e. for which this eigenvalue becomes unstable. The stable orbit is visualised later for the slap 



map 7b (see figure 8(b) i. 

When an eigenvalue crosses —1, a period-doubling bifurcation is expected IIRobi99l 
Chapter 7]. However, in this case no stable period-4 orbit emerges, which would be observed 
as an attractor. We expect that an unstable period-4 orbit is created in this bifurcation; 
presumably it could be rendered stable if the billiard table were modified to have smaller 
curvature around i = 0. 



5.2. Birth of chaotic attractor 

In numerical experiments, a globally attracting chaotic attractor is observed for /I > 0. 1 where 
the period-2 orbit is no longer stable. However, this chaotic attractor is also observed to persist 
for A > A* ~ 0.0712 in coexistence with the attracting period-2 orbit discussed above. 
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Figure 7. (a) Arc length .S4 of the fourth collision point, starting from initial conditions 
xo = (i'o, iPo) at the upper discontinuity sq := — \/3/4 in the cuspless cardioid, as a function of 
the initial angle <po. The dotted horizontal lines show the two discontinuities s = ±\/3/4 and 
4 = (the centre of the vertical line), (b) Basin of attraction of the chaotic attractor for the 
cuspless cardioid with A = 0.09, shown in phase space with coordinates .v and sin((p), where 
(p is the exit angle. The points which are not attracted to the chaotic attractor (shown in blank) 
are attracted to the stable period-2 orbit. The basins appear to be intermingled in some parts of 
phase space. 



How and why does this chaotic attractor arise? For A just above A* it is observed to be 
very thin, and is concentrated around what appears to be a periodic orbit of period 4 joining the 
two points of discontinuity in the curvature (figure |5jb)). We search for such a period-4 orbit 
as follows. Starting exactly at the upper discontinuity, we shoot trajectories with different 
initial angles (p, and follow them for four collisions. Plotting the arc length s of the collision 
point at the fourth collision gives figure 7(a) for angles close to ^ = 0, which are those for 
which the attractor appears. 

We see that for A < A*, the fourth collision occurs for values of arc length greater than s*, 
that is, the collisions occur on the vertical line joining the two discontinuities. The trajectory 
then converges towards the attracting period-2 orbit. For A > A*, the trajectory returns to 
locations with s < s*, that is it collides on the curved surface just above the discontinuity. 
Nonetheless, due to the focusing nature of the boundary, the trajectory is then reinjected into 
the same region, and in this case we observe that the subsequent dynamics converges to a 
chaotic attractor for any A > A*. 

For A = A*, the trajectory returns exactly to the same discontinuity on the fourth 
colhsion, and there is indeed a period-4 orbit joining the two discontinuities exactly at this 
value of A. By searching for the pair (A,(p) for which the curve of figure 7(a) just touches 
s*, we estimate the critical values as A* ~ 0.0712173 and (p* ~ -0.0289796 (see figure [7(a) 1. 
Simulations confirm the existence of the orbit and its instability. 

Since the period-2 and chaotic attractors coexist in the interval A G [A*, Ac], it is of 
interest to examine their respective basins of attraction in phase space, i.e. the sets of initial 
conditions in phase space which converge to each attractor These basins are shown in 
figure 7(b) for the case A = 0.09. In certain regions the basins appear to be intermingled 
arbitrarily closely. 
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Figure 8. (a) Lyapunov exponents for the cuspless cardioid, averaged over 100 independent 
trajectories. Note tliat v+ crosses when the chaotic attractor appears, and has a jump where 
the period-2 orbit becomes unstable, since then there are no longer attracting periodic orbits 
with v+ < 0. (b) Slap billiard map, together with its second and fourth iterates, for the cuspless 
cardioid. 



5.3. Chaotic attractor 

As mentioned above, a chaotic attractor is found for all A > A*, and for all X > Xc this is 
the only attractor. Numerical results again indicate that there is one positive and one negative 
Lyapunov exponent which are both constant almost everywhere, so that the chaotic attractor 
again appears to be non-uniformly hyperbolic. 

As X increases, the attractor expands in phase space, as shown in figure [6] until it fills 
the whole of phase space when X = I. In particular, it visits a larger and larger part of the 
boundary. Note that the "corners" which are visible in plots of the attractors correspond to 
the positions of the discontinuities. Thus the non-differentiability of the map gives rise to 
a non-smooth shape of the attractor The form of the attractor is reminiscent of the Henon 
attractor [Hen76 | and its linear counterpart, the Lozi map IILoz78i . This gives an example 
where the conclusion of IMPS 081 on the existence of complicated dynamics holds, but with 
singularities which are not accounted for in that reference. 



5.4. Slap billiard map for the cuspless cardioid 



Figure 8(b) shows the slap billiard map and its second and fourth iterates for the cuspless 
billiard. The graph of the 4th iterate confirms that there are no period-4 orbits. 



6. Smooth deformations of the circle: three-pointed egg 

The last example of focusing pinball billiards that we consider is a family of smooth 
deformations of the circle, given by the equation 

p(0) = l + acos(30) (13) 

in polar coordinates (p,0) with 9 S [—n^n); they are related to a class of billiards studied 
by Robnik in the case of elastic collisions LRobn83 J . The shape parameter a determines the 
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exact shape of the table: a = gives a circular table, while small values of a give smooth 
perturbations of the circle with 3-fold rotational symmetry, a shape which we call a "three- 
pointed egg", shown in figure 11 (a) However, when a > 1/10 the table becomes non-convex 
and non-focusing, since "waist" regions develop; we thus consider only a G [0, 1/10] in this 
paper. An important feature of this class of tables is that the boundary is analytic, and hence 
contains no singular points. The domain of the pinball billiard map is thus compact. 

As in the ellipse, the arc length is difficult to calculate, and we again replace it by the 
angle 9 of collision with respect to the positive jc-axis. The curvature may be evaluated using 
( 10 1. The minimum curvature occurs for 9o := n/3 and its symmetric counterparts 6i :~ n 
and 02 := — ?r/3; it is 

K{e,)^-\-^^^ for/ = 0,1,2. (14) 

We denote by 0, :— 9i + n (mod 2n) the boundary positions diametrically opposite the 0,; 
these are the points of maximal curvature, with ff(0,) = — ■ 



6.1. Hamiltonian case, X = \ 

It is informative to first study the Hamiltonian case, when A = 1 and the collisions are elastic. 
In this case, for any a £ [0,1/10] the system exhibits a generic mixed-type phase space, 
consisting of elliptic islands and a chaotic sea, as shown in figure [9] for three values of the 
parameter a in the allowed range. 

For a — 0, the dynamics is completely integrable, and the phase space is foliated by 
horizontal invariant curves. As a increases, alternating elliptic and hyperbolic periodic points 
appear ||LL92[ Section 3.2]. Due to the 3-fold symmetry of the tables, there is a bias towards 
period-3 orbits: the figures exhibit elliptic points of period 3, with their associated islands, 
and for small enough a also islands around elliptic period-2 points. For larger values of a, 
there is a transition towards Hamiltonian chaos with a chaotic sea MLL92I Chapter 4], although 
here fuU chaos is never reached. 



6.2. Slap billiard map 

To understand the dynamics of this class of billiards, we next study the slap billiard map 
Tq, which gives information on the periodic orbits which occur when A = 0. Tq is shown in 



figure 10 as a function of the angle around the origin, together with its iterates up to the 
fourth. When the shape parameter a tends to zero, Tq converges to the translation x+n, while 
for a larger than a critical value, three critical points arise, i.e. for which T^{Q) = 0. This 
implies that there are pairs of boundary points whose normals have the same intersection point 
on the opposite side of the table. 

The graph of shows that there are no period-3 orbits when A = 0. Thus the period- 
3 orbits which are found for A = 1 are created at some non-trivial value of A G (0, 1). For 
a = 0.05 there are stable period-2 orbits joining the vertices of the billiard with their opposite 
flat regions. These have become unstable by a = 0.08, and so are no longer visible. There 
are also additional period-2 orbits, which are always unstable - these join parts of the billiard 
at which the tangents are parallel. The following sections study all of these classes of orbits 
in more detail. 
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(a) a = 0.02 





(b) a = 0.05 





(c) a = 0.08 



Figure 9. Phase space of the three-pointed egg in the Hamiltonian case X = I (elastic 
collisions), for three values of the shape parameter a as indicated. Different colours indicate 
trajectories from different initial conditions. The phase spaces are of mixed type, consisting 
of elliptic islands surrounding periodic orbits, together with a chaotic sea in (b) and (c). The 
elliptic period-3 orbits can be clearly seen in all three cases along the lines sin((p) = ±0.5, and 
the hyperbolic ones can be seen in (a) on the same line. In (a) and (b), the elliptic period-2 
orbits can also be observed on the line (p = 0; these have become unstable for the value of a 
in (c). 
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Figure 10. Iterates up to the fourth of the slap billiard map (A = 0) for the three-pointed egg, 
for (a) a = 0.05 and (b) a = 0.08. In (b), critical points have developed, and the previously 
stable period-2 orbits have become unstable in a period-doubling bifurcation, giving rise to 
stable period-4 orbits nearby. There are no period-3 orbits. 
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6.3. Period-2 orbits 



Each vertex 0, of the three-pointed egg is joined to its diametrically-opposite point of 
minimum curvature 0, by a periodic orbit of period 2, with outgoing angle zero; these orbits 



are depicted in figure 11(a) They exist for any a, but can become unstable. Their stability 
matrix, which depends on a and A, is given by 



/ l+a'*-2a^(163+162A,) 
162a^(a^-18a-l)A(A.+l) 



2(a^+18a-l)(l+A) 
(a-l)2 



(15) 



Let jU±(a, A) denote the eigenvalues of the stability matrix. For (a, A) G [0,0.1] x [0, 1] 
we have /i±(a,A) 1. However, defining 



a(A):= 



82 + 162A + 82A2 - 9(1 + A) V83 + 162A + 83A2 

r+A2 



(16) 



we have /i_(ci:(A),A) = —1, i.e. one of the eigenvalues crosses —1, giving rise to a period- 
doubling bifurcation: the period-2 orbits become unstable for a > a (A), and a stable period- 
4 orbit is created nearby in phase space. If a > 6;(0) ~ 0.0781, then the period-2 orbit is 
unstable for all A. The curve a(A) is shown in figure [12] 

On the other hand, for a < a{X), both eigenvalues have modulus less than one, so the 
period-2 orbits are stable. Moreover, for pairs (a, A) such that a < cc{l) ~ 0.0554 and 

,,,, 9(1+A)- V82+160A + 82A2 
a>a(A):=^ ^ — , (17) 



A-1 



the eigenvalues are complex conjugates. 



6.4. Period-2) orbits 

For elastic collisions (A = 1), there are two elliptic period-3 orbits of non-zero angle which 
join the flat parts of the boundary, {0,}, as can be seen in figure [9] These orbits come in 
pairs, with the same points on the boundary but opposite angles, and by symmetry they form 
equilateral triangles in configuration space. It can be verified that they are elliptic (i.e. the 
eigenvalues of the stability matrix lie on the unit circle) for any a. 

When A is decreased away from 1, these period-3 orbits persist, but they become 
attracting, a phenomenon which is well known for weak dissipative perturbations of 
Hamiltonian systems I LT85I IFGH Y96"' 'MG081. Higher-order periodic attractors are also 
expected for very weak dissipation |FGHY96|, and these are indeed also found in our system, 
although we have not investigated this point in detail. 

The symmetry of the table implies that the orbits remain equilateral triangles, but they 
must rotate away from their initial position in order to accommodate the altered reflection 



angles. Figure 1 1(a) shows the change in position of one of these period-3 orbits as A varies 
in the case a = 0.08. 

Below a value A = A^ (a), which depends on a, these orbits cease to exist. To show this 
we use the following argument. Recall that the graph of Tq does not intersect the graph of the 
identity, proving that there are no periodic orbits of period 3 for the slap map. This argument 
extends to for small values of A > 0, as follows. Consider the projection 7i{9,(p) = 9 to 
the space of positions only, and study the set {{9o,noT^{OQ, (p) : (Oq, (p) S M}. This consists 
of all pairs (0o, 63), where 9o is an initial position and 63 is the position of the third iterate 
of the initial point, chosen with any initial angle. If the graph of this set does not intersect 
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Figure 11. (a) Change in position of one of the attracting period-3 orbits as X is varied, 
with the shape parameter a = 0.08. For A = 1 (thick black line) the orbit is elliptic; 
for A < 1 it is attracting. The values of A shown are, in an anti-clockwise direction, 
A = 1.0,0.9,0.8,0.7,0.6,0.5,0.45,0.41,0.40. The latter value is close to the nuinerically- 
determined limit of existence of the period-3 orbits, which is A ~ 0.39. For this value of 
a, the period-2 orbits are unstable; they are shown as dotted lines, (b) Third iterate of the 
perturbed slap map for the three-pointed egg with a = 0.08, for A = 0.0 (slap map; solid back 
line), A = 0.2 (red dots), A = 0.3 (green +), and A = 0.4 (blue x). The fact that there is 
no intersection with the diagonal for A < 0.4 shows that there are no period-3 orbits in those 
cases. 



the graph of the identity map {(0,0)}, then we can conclude that there are no periodic orbits 
of period 3. Of course, if there is such an intersection, then we can conclude nothing - in 
particular, this does not imply that there are such period-3 orbits. 

Figure 1 1(b) shows this map for a — 0.08 and small values of X, from which we indeed 
conclude that Tx has no periodic orbits of period 3 below some A^., with ~ 0.39 in this 
case. For X > Xc, there are intersections, but we reiterate that no firm conclusions can be 
drawn from this; nonetheless, attracting period-3 orbits are found numerically for all X ^ Xc- 

These results suggests that the period-3 orbits disappear in a saddle-node bifurcation at 
X — Xc, by annihilating with another class of period-3 orbits. This latter class of orbits joins 
the vertices 0, of the three-pointed egg when A = 1, and are always unstable. They also rotate 
upon changing A. We conjecture that these orbits collide with the stable period-3 orbits at 
A — Xc- Further evidence for this is obtained by calculating the maximal Lyapunov exponent 
of the stable period-3 orbits as A — > Xc- The Lyapunov exponent converges to 0, suggesting an 
eigenvalue 1 of the stability matrix, which indeed corresponds to a saddle-node bifurcation 
IIRobi99l Chapter 71. 



6.5. Global picture of parameter space 

As we vary the parameters (o:,A) in the two-dimensional parameter space [0,0.1] x [0, 1], we 
numerically observe a variety of different types and combinations of attractors. Roughly, for 
a close to there are only periodic attractors, while for a close to 0. 1 there is a range of A 
for which chaotic attractors are observed to coexist with periodic ones, in particular with the 
period-3 attractors described above. 
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Figure 12. Parameter space of the three-pointed egg, showing the regions of parameters where 
some (red x) or no (green +) chaotic behaviour is found. This was calculated by taking 100 
random initial conditions for each {a, A) pair and calculating the con'esponding Lyapunov 
exponents. If any of the trajectories has a positive Lyapunov exponent, then (a.X) is assigned 
to the region with chaotic behaviour. No chaotic attractors were found for a < 0.06, except 
for the elastic case X = \, where there is always some chaotic behaviour for any a > 0. Any 
attractors whose basins occupy only a very small fraction of phase space are not seen in this 
calculation, and thus do not appear in the figure. The black curve represents the function (x{X) 
calculated in section [63] which gives the limit of stability of the period-2 orbits: above the 
curve, the period-2 orbits are unstable; below it they are stable. 



In order to begin to understand this phenomenology, we first give in figure 12 a global 
picture of parameter space, which depicts only the most gross level of information: whether 
or not chaotic behaviour is found for a given (ct, A) pair, that is whether there are any initial 
conditions which lie in the basin of attraction of some chaotic attractor, characterised by 
having a positive Lyapunov exponent. We can conclude from the figure, coupled with the 
analysis of the stability of period-2 orbits, that increasing a tends to destabilise the dynamics, 
as could be expected. 

We remark that for X — I and any a > there are always some chaotic orbits, namely 
those which wander in the chaotic sea exhibited in figure |9] (or those which stay confined to a 
stochastic layer for smaller values of a, which is the reason why A = 1 is indicated as showing 
chaotic behaviour in figure 12 The red crosses at A = 1 should thus extend over the whole 
range of allowed a. 

Since the billiard table, and hence the pinball billiard map T^, are smooth, we can invoke 
the implicit function theorem to show that periodic orbits, both attracting and hyperbolic, 
should usually be expected to persist under small enough variations both of a and of A . We are 
thus especially interested in the situations where this does not hold, that is when bifurcations 
occur. The persistence of chaotic attractors under perturbation is less clear, since these can 
undergo different types of crisis, where they suddenly change their extent or indeed disappear 
completely due to the interaction with other features of the dynamics IIOtt93l Chapter 8]. 
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Figure 13. Attractors in configuration space for the three-pointed egg with a = 0.08: (a) 
period-4 orbit which persists from the stable period-4 orbit of the slap map (A = 0); (b) period- 
8, after undergoing a single peiiod-doubhng bifurcation; (c): localised chaotic attractor, after 
the accumulation of period doublings; (d): delocalised chaotic attractor, after the merging 
crisis. Here, the trajectories tend to remain for a long time in each part of the attractor that was 
previously localised, before jumping to a different part, as shown by the colours in the figure. 
In each of (c) and (d), a coexisting period-3 orbit is also shown in black. Note that for this 
value of a, the horizontal period-2 orbit is unstable even for A = 0, and thus does not attract 
generic initial conditions. 



6.6. Period-doubling cascade 



The main feature visible in figure 12 is the "tongue" region for large a where chaotic attractors 
occur. To understand this, we consider a cut across parameter space by fixing a = 0.08 
and studying the evolution of the numerically-observed attractors as A is varied. This gives 
representative results for other values of a in this region. 

As discussed above, for this value of a, there are attracting period-4 orbits for A = 0, 
and numerically these are found also for small A. For A > 0.39, there are attracting period- 



3 orbits, as shown above. But figure 12 shows that there is also a region of values of A, 
approximately the interval [0.42,0.66], which overlaps the region of existence of attracting 
period-3 orbits, where chaotic attractors are also found. It is thus of interest to investigate the 
route by which the chaotic behaviour is first created and then destroyed. 

Numerically we find that this chaotic behaviour arises by the well-known period- 
doubling route IIOtt93l Chapter 8], where a cascade of an infinite number of successive period 
doublings at increasing values of A accumulates at a certain value Ax,. For A > Aoo, there 
is a chaotic attractor with a positive Lyapunov exponent, except in windows where periodic 
behaviour is observed, similar to the periodic windows observed in the logistic map iOtt93l 
Chapter 2]. The development of the observed attractors is shown in configuration space in 



figure 13 and in phase space in figure 14 



For small values of A, there are only stable period-4 orbits. At a critical value of A, 
these become unstable via a period-doubling bifurcation, giving stable period-8 orbits. As A 
increases, we observe a repeated doubling of the period of the attracting orbits, giving orbits 
with period 2" at A„, where A3 ~ 0.372, A4 ~ 0.389, A5 ~ 0.393 and Ae ~ 0.3935. These 
values are slightly sensitive to the tolerance used in the Newton-Raphson algorithm, and to 
the precision of the calculation. To obtain these orbits numerically, it is useful to note that 
they can be distinguished from the coexisting period-3 orbits by the respective values of their 
maximal Lyapunov exponent. Although this is always negative for periodic orbits, its value is 
around —0.3 for the attracting period-3 orbits, and around —0.01 for period-32, for example. 

By A = 0.3937 a chaotic attractor (i.e. one with a positive Lyapunov exponent) is 
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Figure 14. (a) Attractors in phase space for the three-pointed egg with a = 0.08 and X = 0.43 
(red), X = 0.45 (green). The coexisting period-3 attractors are also shown (black x). (b) 
Enlargement of the attractor for X =0.45. 



X 




Figure 15. Part of the bifurcation diagram of the three-pointed egg with a = 0.08, for a range 
of X including the periodic-doubling cascade, and the subsequent crisis of the chaotic attractors 
when three spatially distinct attractors combine into a single one. For each horizontal slice 
with fixed X, attractors in phase space obtained from different initial conditions are shown, 
indicating the range of possible behaviour for each value of X . The "noise" visible for smaller 
values of X are period-doubling cascades of other orbits, which also reach a crisis when they 
are absorbed into other attractors. 



obtained. This chaotic attractor has features which appear visually similar to the type of 
structure found in the Henon attractor IIHen76ll . Part of the bifurcation diagram is shown in 
figure[T5]for a — 0.08, where the period-doubling cascade and the subsequent merging crisis 
(see below) can be seen. 

Also visible on figure 12 are what appear to be isolated regions with chaotic attractors. 
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The fact that they appear isolated is apparently due to the resolution of the plot. These chaotic 
attractors also appear to be created by period-doubling cascades, starting from periodic orbits 
of higher period, as is visible as what initially appears to be "noise" in figure 15 In this work 
we have not attempted to exhaustively characterise the structure of attractors, which are known 
to exhibit very complicated coexistence phenomena in other models IIFGHY96I IMG08I . 



6.7. Crises of chaotic attractor 

We continue to fix a = 0.08. Above Ac — 0.42, the attractor has behaviour reminiscent of 
the logistic map, with alternating chaotic and periodic windows. As A increases, the chaotic 
attractor expands in phase space, and has an appearance somewhat similar to that of the Henon 
attractor 

Initially there are three separate attractors, each concentrated near one of the diametric 
period-2 orbits. However, at around A = 0.4435, there is a first crisis ||Ott93l Chapter 8] in the 
chaotic attractor: the previously separate three symmetric parts of the attractor abruptly join 
together into a single attractor. 

This single attractor then continues to expand in phase space until another crisis occurs: 
the chaotic attractor abruptly vanishes when A ~ 0.66. The reason for this is found in 



figure 16(a) which shows the basins of attraction of the chaotic and period-3 attractors for 
A — 0.66, just below the crisis. As A increases, both the chaotic attractor and the basins of 
the periodic attractors grow. For a slightly larger value of A, the chaotic attractor touches the 
basins of the period-3 attractors. After a chaotic transient which resembles the previously- 
existing attractor, all initial conditions are then attracted to the period-3 orbits, and the chaotic 
attractor thus disappears. 



In figure 16(a) the basins of the two distinct period-3 orbits are not distinguished. In fact, 
they are completely separate: the points close to each orbit are attracted to it, and the basins 
are divided by that of the chaotic attractor. Once the chaotic attractor has vanished, however, 
the region that was occupied by its basin is now attracted to one of the two period-3 orbits. 



Their individual basins are shown in figure 16(b) We find that their basins are apparently 



completely intermingled in this intermediate region. 

We remark that the chaotic attractor found in the three-pointed egg is interesting since 
it is the only case that we have studied of a billiard table having a smooth boundary (that is, 
with no points at which the boundary is not C^) which exhibits a chaotic attractor. 



7. Discussion and conclusions 



We have shown that there is a large variety of interesting phenomena in dissipative pinball 
billiards. Generically, a mixture of periodic and chaotic attractors is found, with only the 
former in the ellipse and only the latter in the cardioid. 

In the cuspless cardioid, we showed that discontinuities in the curvature of the table 
boundary can lead to the the sudden creation of a chaotic attractor These discontinuities are 
also reflected in the existence of a kink structure in the resulting attractor in phase space. 

Smooth tables, such as the three-pointed egg, can exhibit complicated sequences of 
bifurcations, including period-doubling cascades from periodic orbits of different periods. 
These can lead to chaotic attractors which then collapse by means of crises when they touch 
a different basin of attraction. 

In MMPS0 8 1, it was shown that in systems with dominated splitting there is another 
possibility in addition to periodic and chaotic attractors, namely the existence of sets 
diffeomorphic to a circle where the dynamics is conjugate to an irrational rotation (for details. 
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(a) A = 0.66 (b) X = 0.7 

Figure 16. Basins of attraction in phase space for the three-pointed egg with a = 0.08 and (a) 
A = 0.66; (b) A = 0.7. In (a), the chaotic attractor is shown with blue points and its basin of 
attraction with green +; the period-3 attracting orbits are shown as black x , and their basins of 
attraction in white. The individual basins of the two period-3 orbits are disjoint, with one lying 
above and the other below the chaotic attractor. In (b), the chaotic attractor has disappeared 
in a crisis, leaving just two attracting peiiod-3 orbits, shown as black x . The basins of each 
of these orbits are distinguished, one in red dots and the other blank. The region that in (a) 
was occupied by the basin of the chaotic attractor is now a region where the basins of the two 
periodic orbits intermingle, apparently arbitrarily closely in the whole region. 

see IIMPS08I Section 2]). We have not observed such behaviour in the models we have studied. 
One of the possible reasons for the absence of such behaviour is the fact that we have only 
used uniform contraction of the angle at collisions; other types of collision rule are under 
study. This kind of behaviour could, however, perhaps be found on the boundaries which 
divide different basins of attraction, which are often found to be smooth one-dimensional 
curves, as depicted in several of the figures in this paper. 

We found that weakly-dissipative pinball billiards, i.e. with A close to 1, give some 
information on the existence of elliptic periodic orbits in the elastic case. Furthermore, we 
studied two cases, the cardioid and the cuspless cardioid, for which the dynamics is fully 
chaotic when A = 1, and in this case weak dissipation leads instead to a chaotic attractor. We 
suggest that adding weak dissipation could thus possibly be useful as a means of detecting 
small elliptic islands in the phase space of Hamiltonian billiards. An interesting alternative 
method for doing this was proposed in ref. IITK07I . 

An open question is whether a similar dynamical portrait in parameter space occurs 
for similar billiards, for example a /T-pointed egg given by the polar equation p{9) — 
1 + acos{KO), for other values of K > 3 eN. 

The study of pinball billiard dynamics in other types of billiard table, including 
dispersing billiards, is a work in progress. 
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Appendix A. Numerical methods 

In this appendix we summarise the numerical methods used to calculate the figures shown in 
the paper. 

Appendix A.l. Trajectories and attractors 

Initial conditions for the simulations are generally chosen uniformly in the interior of the 
billiard table, with initial velocities chosen with uniform directions <p e [0,2;r), and are then 
mapped to the first intersection with the table boundary. The subsequent dynamics is followed 
by finding each consecutive intersection with the boundary and implementing the pinball-type 
colUsions. Cartesian coordinates are used whenever possible. 

The distance (or equivalently time) f until the next intersection with the table boundary 
is found using techniques specific to each type of table. For the ellipse, the intersection is 
given by solving a quadratic equation which is a generalisation of that used for intersections 
with a circle IIGas98J . In the cardioid, a Newton-Raphson step is used to find the first 
boundary collision, and thereafter a cubic equation is solved, following |BD97|. In the 
cuspless cardioid, this is adapted so that if the proposed root lies in the forbidden region 
then the particle instead collides with the plane. In the three-pointed egg, a Newton-Raphson 
method is used to find the intersections. 

Once the collision point q on the boundary has been found, the pinball billiard dynamics 
is implemented as follows. The normal vector n(q) at the point q of collision is calculated 
based on the equation describing the table boundary, and a tangent vector t(q) there is 
constructed orthogonal to n. The incoming angle r\ at the collision is found as Tj = 
cos^'(v-n), and the outgoing angle is then given by (jO = Xr\. The outgoing velocity in 
Cartesian coordinates is finally given by v' = ncos (jO + tJsin^, where d — ±\ ensures that 
the component of the outgoing velocity is in the same direction with respect to the tangent 
vector as the incoming velocity was, as determined by the sign of v • t. Note that in the case 
of elastic collisions, the dynamics may be implemented completely in Cartesian coordinates 
|[Gas98 |; however, the definition of pinball billiard dynamics makes the use of trigonometric 
functions unavoidable. 

To find attractors of the dynamics, the first several thousand collisions in the trajectory 
emanating from each initial collision are discarded. Assuming that transients have already 
decayed by that point, and that the dynamics is transitive on the attractor, the result is an 
approximate picture of the attractor. 

Appendix A.l. Lyapunov exponents 

The Lyapunov exponents of an orbit of the pinball billiard map Tx are calculated using the 
algorithm of Dellago et al. [DPH96|. This involves the derivative Q of the map of the pinball 
billiard transformation. For completeness we summarise the main steps of the procedure. 

Given an initial point x in the two-dimensional phase space, several thousand collisions 
in its trajectory are discarded to remove transients. Then the trajectory is extended to the next 
collision, and two orthonormal two-dimensional vectors vi and V2, which are thought of as 
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tangent vectors at the point x, are manipulated in parallel. Each tangent vector is acted on by 
the derivative matrix DJx , to give two new vectors v j and Vj, which are then reorthogonalised 
using a Gram-Schmidt procedure. This process is repeated at each collision. The Lyapunov 
exponents Vi and V2 are given by the average growth rates of the normalising factors. In 
particular, their sum Vi + V2 gives the logarithm of the mean rate of volume contraction of the 
dynamics on the attractor in phase space. 

Appendix A.3. Basins of attraction 

To determine basins of attraction in phase space, initial conditions are chosen uniformly 
in phase space as above (with the initial boundary collision being an elastic one, to allow 
all initial directions to be produced). For each initial condition {s,<p), the two Lyapunov 
exponents are calculated. If the maximal Lyapunov exponent is positive, then the trajectory 
converges to a chaotic attractor; if it is negative or then it converges to a periodic 
attractor. A symbol is plotted at {s, (p) according to the type of attractor found there. For 
the correspondence between negative Lyapunov exponents and periodic orbits, and between 
positive Lyapunov exponents and chaotic behaviour on strange attractors, we refer to IIER85I 
Section III]. 
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